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Abstract: 

We perform a calculation in one-loop chiral perturbation theory of the two-pion matrix elements and 
correlation functions of an / = 0 scalar operator, in finite and infinite volumes for both full and quenched 
QCD. We show that major difficulties arise in the quenched theory due to the lack of unitarity. Similar 
problems are expected for quenched lattice calculations of AT —> tttt amplitudes with AI = 1/2. Our 
results raise the important question of whether it is consistent to study K —>■ tttt amplitudes beyond 
leading order in chiral perturbation theory in quenched or partially quenched QCD. 


PACS numbers: 11.15.Ha,12.38.Gc,12.15Ff 


1 Introduction 


A precise quantitative evaluation of weak non-leptonic amplitudes in kaon decays is an enormous challenge 
for lattice QCD. Although it has been demonstrated that such a calculation is possible in principle, 
a number of major practical difficulties must first be overcome. These difficulties are related to the 
construction of finite matrix elements of renormalized operators from the lattice bare ones and to the 
extraction of physical amplitudes, including final state interaction phases, from Euclidean correlation 
functions. For the latter problem, it has been demonstrated that it would be possible in principle to 
obtain the physical amplitudes by performing unquenched simulations with physical quark masses on 
lattice volumes large enough to have discretization errors and finite size effects under control [1, 2]. 
At present however, it is not possible to perform unquenched simulations on such large volumes and 
therefore a certain number of approximations are necessary. One of the main approximations (in addition 
to quenching) consists in working with unphysical quark masses and/or external meson momenta, and 
estimating the physical amplitudes by extrapolating to the physical point. A key element of our strategy 
in evaluating K tttt matrix elements is the use of Chiral Perturbation Theory (yPT) at next-to-leading 
order (NLO) [3]. In a recent paper [4], we have presented the relevant formulae for A/ = 3/2 transitions 
on finite and infinite volumes, in the full theory and in the quenched approximation. Our results show 
explicitly that all corrections which vanish as inverse powers of the volume can be eliminated by using 
the methods introduced in refs. [1, 2] The remaining finite volume corrections are exponentially small 
(of At the NLO in yPT this is true also in the quenched approximation. 

In this paper we present the results obtained at NLO in yPT for matrix elements with an / = 0 
two-pion final state. We study these matrix elements in order to illustrate the main features present 
in AI = 1/2 K TTTT transitions in quenched QCD. For these decays, the lack of unitarity of the 
quenched theory leads to a number of problems which need to be solved in order to understand the 
volume dependence and to extract the amplitudes. The main consequences of quenching, due to the lack 
of unitarity, can be summarized as follows: 

• the final state interaction phase is not universal, since it depends on the operator used to create 
the two-pion state. This is not surprising, since the basis of Watson’s theorem is unitarity; 

• the Liischer quantization condition [6] for the two-pion energy levels in a finite volume does not 
hold [8]; 

• a related consequence is that the Lellouch-Liischer (LL) relation between the absolute value of the 
physical amplitudes and the finite volume matrix elements [1, 2] is no longer valid. It is therefore 
not possible to take the infinite volume limit at constant physics, namely with a fixed value of W ; 

• whereas it is normally possible to extract the lattice amplitudes by constructing suitable time- 
independent ratios of correlation functions, this procedure fails in the quenched theory as explained 
in section 2.2. In particular the time dependence of correlation functions corresponding to different 
operators which create the same external state is not the same; 

• in addition to the usual exponential dependence on the time intervals, the presence of the double 
pole corresponding to the incomplete r]' propagator generates, at one-loop order in yPT, terms 
in the Euclidean correlation functions which depend linearly, quadratically or cubically on the 
time [7]. Unlike the corrections which shift the two-pion energy in a finite volume [4], these terms 
do not exponentiate and may cause practical problems in the extraction of the finite volume matrix 
elements. A related problem is the appearance, at fixed L, of corrections linear or cubic in L [7, 8]. 

The dependence on the r]' parameters can be removed by working in partially quenched QCD, where 
the ry' is heavy and decouples from the light Goldstone boson sector, although, in general, residual double 
poles will remain. All the problems originating from the lack of unitarity (denoted as the unitarity problem 
in the following) would, however, remain the same. Unitarity is recovered from partially quenched QCD 

^The infinite-volume limit, L ^ oo (where L is the length of each spatial dimension of the lattice), is to be taken at 
fixed physics, i.e. at fixed two-pion energy W. 
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only in the limit when the number of sea and valence quarks is equal and their masses are equal. This 
corresponds to full QCD. 

In view of the difficulties listed above, it may be questioned whether it is possible to obtain K —>■ tttt 
decay amplitudes beyond leading order in yPT in quenched or partially quenched QCD. In the absence of 
a solution to the problems encountered and discussed in this paper, we would be limited to extracting the 
effective couplings (the low-energy constants) corresponding to the operators in the weak Hamiltonian 
at lowest order in the chiral expansion from K ^ n matrix elements computed in lattice simulations. 
The absence of unitarity is intrinsic to quenched and partially quenched QCD, and we do not have a 
solution to the unitarity problem. Nevertheless it is tempting to speculate whether there might not be 
a possible pragmatic way to proceed, in spite of the failure of Watson’s theorem. Indeed the key point 
is that, in the quenched case, the two-pion state is no longer an eigenstate of the strong interaction 
Hamiltonian. The eigenstates of the Hamiltonian would be, formally, linear combinations of physical 
pions and unphysical mesons composed of the pseudo-fermion fields. The latter, however, have the wrong 
spin-statistic properties and for this reason unitarity breaks down. We speculate that it might be possible 
to recover a variation of Watson’s theorem, and of the LL formula in finite volumes, by a suitable re¬ 
interpretation of the quenched theory, for example by using the replica method of ref. [5], working in the 
basis of Hamiltonian eigenstates. However, we stress that this is only a speculation and we will report 
on the conclusions of our investigations of this important question in a future paper. 

Since all the major difficulties arising from the unitarity problem depend only on the quantum numbers 
of the operators and final state, for the sake of illustration we discuss in this paper matrix elements of the 
form (7r7r|S'|0), where S' is a scalar and isoscalar operator which can annihilate the two-pion state. We also 
discuss the properties of the correlation functions from which such matrix elements are obtained. The 
discussion can readily be extended to K ^ mr matrix elements of AI =1/2 operators of the effective weak 
Hamiltonian [9]. The results for the AI = 1/2 amplitudes in one loop yPT on finite and infinite volumes, 
in the full theory and in the quenched approximation, will be presented in a forthcoming publication [10]. 
In our calculation we have used the formulation of the quenched chiral Lagrangian introduced in ref. [11], 
using the conventions and notation presented in section 3 of ref. [4]. The scalar operator is defined by 


S' = tr 




( 1 ) 


in the full theory and as 


S« = 


1 = 1,3 


t 


a 


( 2 ) 


in the quenched theory, where the trace (sum) is taken over the indices of the chiral group S17(3 )l 0 
SU{3)r (graded group SU{3\3)l <8> S17(3|3)fl). In the quenched case the field S'^ is the graded extension 
of the standard field S of the full theory. 

The main results of our calculations are presented in the four appendices A~D, which contain the 
following: 


1. the NLO expression for the matrix element of the scalar operator in full QCD and in infinite volume, 
(7r+7r“|S|0). This is given in eq. (14); 

2. the expression for the corresponding quantity in the quenched theory (7r+7r“|S*|0), presented in 
eq. (19); 

3. the NLO result for the correlation function (0|7r/j-(ti)7r/r(t2)S(0)|0) in a finite volume and in full 
QCD. This is given in eq. (20); 

4. the corresponding correlation function in the quenched theory, (0|7r/j-(ti)7r/r(t2)S"^(0)|0), presented 
in eq. (26). 

In the above S'(O) = S'(ai = 0, t = 0) with the corresponding definition of 5''*(0). The expressions for the 
matrix elements are given in Minkowski space, whereas those for the correlation functions are presented 


2 



in Euclidean space. In the correlation functions we have used the following definition for the Fourier 
transform of the fields 



d^x Tr{x, t) . 


( 3 ) 


The correlation functions depend on the choice of ti and In this paper, for purposes of illustration, 
we present the results for the two cases ti = t 2 and ti ^ t 2 - 


In section 2 we discuss the physical interpretation of the expressions obtained at one-loop in yPT 
for the different cases and the implications for the relation between finite volume Euclidean correlation 
functions and physical amplitudes (including final state interaction phases) [1, 2]. The relevant expressions 
and the technical details can be found in the appendices. 


2 Discussion of the One-loop Calculations 

In this section, we discuss the extraction of the physical amplitudes from finite-volume Euclidean correla¬ 
tion functions, using the results obtained in one-loop yPT- We first consider the unquenched case, where 
we explicitly check the validity of the LL relation [1, 2], derived using general properties of quantum 
mechanics and field theory. 

2.1 Extraction of the Physical Amplitude from the Scalar Correlation Func¬ 
tion in Pull QCD 

We begin our discussion from the correlation function of the scalar operator with two pion fields, in the 
full theory at finite volume, (0|7r^-(ti)7rr(t2)S'(0)|0). The tree-level and one-loop diagrams are shown in 
figs. It and Iz, la and lb. Final state interactions, and consequently power-like finite volume corrections, 
are only given by the diagram in fig. lb The NLO expression for the correlation function is given in 
eq. (20). The corresponding Minkowski amplitude in infinite volume is given in eq. (14). In this case 
we denote the contribution of the diagrams in fig. Iz, la and lb by 1^, la and Ib respectively and define 
the relative one-loop correction to the infinite volume amplitude as Aod = l/(47r/)^ {Iz + la + h), see 
eqs. (14) and (15). Final state interactions are encoded in the function A{m) introduced in eq. (16). 

In a finite volume the corrections to the amplitude from the diagrams in figs. Iz and la {Iz and la 
respectively) are the same as in infinite volume up to exponentially small terms (in the volume) that will 
be neglected in the following. The diagram in fig. lb gives a correction, Ib{ti,t 2 ), which is a function of 
the time coordinates of the interpolating operators which annihilate the two pions (the two-pion sink). 

At lowest order (obtained by setting Iz = Ia= /&(ti, ^ 2 ) = 0), the time dependent factor 

- -^, ( 4 ) 

can be removed by dividing by the two-pion propagator in the free theory, and in this way the required 
matrix element can be obtained. In equation (4) E is the energy of each of the pions, E = 

At one-loop order the relevant finite-volume correction to the correlation function of eq. (20) is therefore 
given by Ib{ti,t 2 ) which is presented explicitly in eq. (21). We rewrite Ib{ti,t 2 ) as Ib{ti,t 2 ) = I4e(/t,) -I- 
T{t 2 ) +IZ{ti,t 2 )t where Ib is the corresponding infinite-volume one-loop contribution to the matrix 
element (see eq. (15)) and we now discuss the significance of the terms T and IZ. 

T{t 2 ) contains the one-loop corrections which are multiplied by the correct time dependence, 
exp(—lTt 2 ) exp(—i?(ti — ^ 2 )) (after exponentiation), where W = 2E AW and AW is the shift of 
the two-pion energy due to interactions in the finite volume [4]. We can readily extract AW from the 

^The evaluation of this diagram for 1 = 2 final states was explained in some detail in section 4 of ref. [4]. We therefore 
do not present a description of the calculation for 7 = 0 final states, but limit the discussion to the implications of the 
results. 
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Figure 1: Tree level (t) and one-loop xPT' diagrams for the (tt 7r+|S'|0) amplitude and (0|7r 7r+S'|0) 
correlation function. The grey circle represents the scalar source while the sguares are strong vertices. 


coefficient of t 2 in the expression for T{t 2 ) in eq. (23), 




(5) 


It is then straightforward to check that eq. (5) reproduces the well known result for the scattering length 
[13]. When the two pions are at rest eq. (5) gives 


W = 2to^ 


47ra^=° 


( 6 ) 


where Oo“° = 7m^/(167r/^). 


T(t 2 ) also contains the finite volume corrections to: 

1. the matrix element of the scalar operator. These are given by the one-loop component of the 
LL-factor relating the infinite-volume and finite-volume amplitudes, AIoo and My respectively, 

= LL X IMvH; 


2. the two-pion sink used to annihilate the pions created by the scalar source. We refer to these as 
Forward Time Contributions or FTCs. They are presented explicitly in eq. (25). 


The FTCs are eliminated by dividing (0|7r^-(<i)7r-(t2)*S'(0)|0) by (the square root of) a suitable tttt 
correlation function [4] 

G777r^777r(tl,t2) = ^ (0|7r^(tl)7r_^(t2)7r])-(-t2)7rt(-ti) |0) (7) 

p.f; |p|=l9l. 171 fixed 


o-2lVt2 


{2EY 

^—2Wt2—2E{t\—t2) 


1 + 1 - 


3m^ \ 2Z. \ 

8F;2 j 3pEL^ ) 

3ml 


L° 1+1- 




8E^ J 3pEL^ 


for ti = t 2 


for > ^2 • 
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Note that in the unquenched case the finite-volume energy W appearing in eq. (7) is, as expected, the 
same as in 

(0|<^(ti)^-(t2)^(0)|0) a , (8) 

because we are considering the same final state. As shown below, this is not true in the quenched and in 
the partially quenched theories. 

Finally TZ{ti,t 2 ) corresponds to contributions to the correlation function whose time dependence is 
governed by energies different from W and which therefore can be eliminated by studying this time 
dependence [4, 12]. Since it is not relevant to our discussion, TZ{ti,t 2 ) will be neglected in the formulae 
given in the appendices. 


2.2 The Scalar Correlation Function in the Quenched Theory 


In one-loop yPT with 1 = 2 final states (as for example in A/ = 3/2 AT ^ tttt transitions) we can 
follow the path outlined above for full QCD also in the quenched theory. This is not the case however, 
when we consider the infinite-volume amplitude or the finite-volume correlation function of the isosinglet 
scalar operator with two pion fields in the quenched theory. The one-loop expression for the amplitude 
(tt+tt" I S"? I 0) is given in eq. (19) and for the correlation function (0 17r/^-(fi)7rZ(t2)5'*(0) | 0) in eq. (26). 
In this case, even in infinite volume there are severe difficulties, for example the imaginary part of the 
amplitude given in eq. (19), corresponding to the final state interactions, depends on the scalar operator 
and diverges at threshold, i.e. when s —> 4m^ [14, 15]. We now consider what happens for the correlation 
function in a finite volume. 

For illustration, let us start by neglecting the contribution of the ij' double pole by setting toq = 0 
and a = 0 (they are defined in eq. (18)). For simplicity we also neglect the terms proportional to vi ^2 
(vi .2 are also defined in eq. (18)). This is similar to the situation encountered in the partially quenched 
case where the rj' is heavy, but unitarity is violated. In this case all the a'll defined in Appendix D satisfy 
a^l = 0 and for E < rriK < rriss only the first line of T‘^{t 2 ) in eq. (29) contributes 


r«(t2) 


E^ 

\ + 1 


w 

f;2l3 ^ ' 

1 3A3L3’ QE^L^J 


^(0) ^(1) Y 

2E^L^ 2tt'^El) 


(9) 


The terms (...;...) are the forward-time contributions {ETCs) respectively for the two cases ti ^ t 2 and 
t2 =ti. 

The first term in eq. (9) is the shift in the energy of the two-pion final state. This shift. 


AIF = - 


/2L3 


is different from that obtained from the tttt ^ tttt correlation function in eq. (7), which is 


AW = - 


V 

ApL^ 



( 10 ) 


( 11 ) 


The latter in fact is the same in the quenched and unquenched theories, since, at this order, it is given 
by a tree diagram. Once the contribution from the r]' double pole is included, at this order there are 
quadratic and cubic terms in O present in T‘^{t 2 ), but not in the tttt ^ tttt correlation function. 

The second term in eq. (9) should be cancelled when extracting the matrix element of the scalar 
operator by dividing by the square root of the tttt ^ tttt correlation function: 


|(7r7r|S"?|0)y| 


'k+-(^i)7r,-(O)S"^(0)|0) 




TTTT—^TTTT 


(G, 0) 


( 12 ) 


since it is the finite volume correction to the sink operator used to annihilate the two pions. One can 
readily verify that the cancellation does not occur, unlike in the unquenched theory. Thus the power 
corrections in 1/L to |(7r7r|S'‘^|0)y| are not those expected on the basis of the Lellouch-Liischer formula. 
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Conclusion 


In this paper we have shown that the standard strategy for extracting the amplitude and the relative 
phase from finite-volume calculations of correlation functions fails in the quenched theory for 1 = 0 two- 
pion final states. We encounter the same problems when computing the matrix elements of the operators 
relevant for AJ = lj2 K ^ tttt transitions at one-loop order in yPT [10]. 

At present we do not know whether it might be possible to recover a modification of Watson’s theorem 
and of the LL formula in finite volumes, which would allow for a consistent determination of / = 0 two-pion 
matrix elements beyond leading order in quenched yPT. We are currently investigating the possibility 
of performing a suitable analytic continuation in the numbers of flavours, using the replica method of 
ref. [5] at one-loop order in yPT. 
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Appendices 


A The Scalar Amplitude in Full QCD 

To fix the conventions we start by writing the chiral Lagrangian used in our calculations in the full QGD: 


/:strong = y tr [(a^Et) (a^E) + Etx + X^E ] , 

where the decay constant / is normalized in such way that ~ 132 MeV. 

The one-loop matrix element of the scalar operator in full QGD in infinite volume is given by 

1 


(13) 


where 


(tt (( f ) 7 r +(- g ) jAjO ) =-^2 [1 + - 4 oo ] =-^2 


1-k 


(47r/)' 


ih +Ia+ h) 


h = - 


|„4 log + I (ml, + 2m|) log , 


m m \ 5m: 


(14) 


h = {ml - 2s)A{m^) - sA{mK) - ^A{mn) 


+ ( ^ml - 2s 


)log(^). log (A) 


4 - ^ 


3s — —mt 


(15) 
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In the above expressions /i is the renormalization scale, s = (p^- +p^+)^ is the square of the two-pion 
center of mass energy and 


A{m) 



(16) 


To the one-loop corrections listed in the above equations, we should add those proportional to the 
Gasser-Leutwyler coefficients, [13], appearing at NLO in the strong interaction chiral Lagrangian. 
These effects, combined with Iz, are reabsorbed into the renormalization of the decay constant, leading 
to the replacement of the factor 1//^ by 1 /in the matrix elements of the scalar operator. Since these 
terms do not affect the finite volume corrections which are the object of the present study, we will not 
discuss them further here. Similarly the matrix elements of the (higher-dimensional) counter-terms of 
the scalar density do not contribute to the finite volume effects and will not be considered. 

The s-wave phase shift for the 1 = 0 two-pion state at one-loop in yPT is given by 


d(s) 


2s — 
167r/^ 



(17) 


For s > 4m^, we have Arg [(tt ((f)7r+(—(f)|S'|0)] = (l(s). Note that the amplitude remains finite at 
threshold, i.e. as s ^ 4 to^. 


B The Scalar Amplitude in Quenched QCD 


In this appendix we discuss the matrix element of the scalar density in quenched QCD. The quenched 
chiral Lagrangian used in our calculations is: 

-mo$o-h a(9p«'o)(9^d>o), (18) 

with the super-77' field, <I)o = f/2i str[logE^]/-\/6. 

The one-loop matrix element of the scalar operator in the quenched theory, at infinite volume, is given 


by 

where 

/« 

Z 

n 


n 


^+^-|59|0)^-_[1 + A^J = -_ 


1-f 


(47r/)' 


{II+ 11 +ID 


= 0 , 


3 V M 
3 


m 


K 


log 


= s 


-f 


1 


- A{mD - ^A{mK) ) + ( ^ml - s ) log 


2 ' "' 2 
2m1 {V2 - 77i) -I- SUi 


2 


2mQ 


1 + 


log ( 


mi: 


ml: 


\ 


+ 


1-f 2 


1 


-m 


+TOn 


3 

4m^ 


A{mss) + 2A{mD + 2 log 


Cf) 


3(s — 4 to2 ) 


A{mD 


4m:: 


log ( 




- 1-f 


3 

-flog 

s — 5mJ 
s — Iml 


K 


9 

mi. 


A{mz 


+ml 


(Aio - A { mss )) - 


8m: 


18 (m^ — c/oyiii^j^ 

2 rQ0..^6 I 1 79^2 _ „'i _ s.™2 „2 i „3 qo^2 ^2 79.^2 


9s(s — 4m^) 9s{m‘D — rn^) 


log 
- 3 

log 


m' 


K 


m^[32m® + 16m|j-(2m^ ~ s) ~ 8m^s^ -f s^ — 32m|-mJ(2m^ — s)] 


18s{m‘D — m2)2(s — 4m^)^ 


-A(m^) 


(19) 


log ( 


m: 


2 \ 1 
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with 




16m^ 


9s(s — 4m2) 9(m|. — )2 

-.2 ',2 


[rnl^A{mss) - 


'm^{2m\ — s) ( m 


9s(m|. — m^) 


mJ[8(TO|. — to^)^(4to^ + 2mjs — s^) + m^s{s — 4to^)^] 


+ - 


2mi 


1 - 


9s(m|. — m^)^(s — 4m2)^ 
4mi TO? 


21(to^) 


3s(s — 4 to2 ) 12(to|. — to^)2 


(mlAio - ml^A{mgs)j 


[8(to|- — to^)^(4to® + 22to^s — IOto^s^ + s^) + m'^s{s — 4to^)^] 
12s(to|- — to2)2(s — 4to2)2 


m’^rrigg + s(4to|- — 5to^) 
Qs{m\ — TO^) 


x.sfe'l- 


V A 


A{s-ml) 


fi i 2 2^1og(^ 


log 


A(to^) 


A 


= 2to|- - TO^ , 


10 


= A log 


l_2l^+A 

l-2ll^-A, 


- m9 (Vs - {mgs + TO^)) , 


m 


A = Vl + 4^+4^-4^(l + 2^). 

S'^ S'^ s s 


The imaginary part of the amplitude diverges at threshold, i.e. as s ^ 4 to^ [14, 15]. 


C Finite-Volume Scalar Correlation Function in Full QCD 


In this appendix we give the complete one-loop expression of the finite-volume correlation function of the 
scalar operator with two pion fields in the full theory: 


(0|7r+Vti)7r-(t2)S'(0)|0) 




1 

(47r/)2 


(-^2 +-^a) H"^ 2 ) , 


( 20 ) 


where E = \J ^ A to^ . Since we are interested in finite volume corrections we give only the explicit 
expression for /&(<!, 12 )- We write 


where 


^00 


All 


A 22 


h{ti,t2) 


W 


[Aoo -I- All + A 22 ] , 


( 21 ) 


1 


1 


6£;2 ^ wl \ 2{E - Wo) 2{E + wq) 

k 


—e 


2{E—wo)t2 


1 


1 


2{E — Wo) 2wo 

J_ f d-r(wi) _ d-(wi) 

L3^i2(£;-wi) 2(£;-twi) 

k 


^2Et2 — 2wot 


1 fJ: - 1 - 

\2wo 2{wo + E) 


_ 2(E-w,)t2 ( d+jwi) do{wi) \ 2Et2-2w,U ( Mwi) _ d-{wi) \ 1 

\2{E — wi) 2wi ) \ 2'u;i 2(wi-I-if)/ j 

J_ f C+(W2) _ C_(W2) 

T3 2^ \2(£;-W2) 2{E + W2) 

k 

_ 2(E-W2)t2 ( C+(W2) Co(w2) \ 2Et2-2w2t^ ( ^ 0 (^ 2 ) _ C-(w2) \ 1 

\2{E — W 2 ) 2w2 ) \ 2w2 2{w2-\-E))1 



In the above expressions 


Wo 


= \ k'^ 


SS 5 

2E'^±Ew + w^ 2 1 j ^ o / \ 

c±H = 3- -^2 -; coM = -^; d±,oH=2c±,oM-2;^. 


wi = yfc^ + ml ; W2 = y k"^ + 


'-K > 


Evaluating the sums, we find 

(0|7r;-(<i)7rr(t2)S'(0)|0) = 
where, for E < mx < rngs, 


^ — Eti ^ — Et2 
“ 2 ^ 


t^2 

nt2) = 


Vt2 


2/2 [ E2L3 


”72 ) R-e(“ 2 loo) + '^(^ 2 )] , 




3E3L3’ 6E3L3, 


/ 2 :( 0 ) _ 3w^ , 

\ 77-3 7- 3 V O 7712 / 


2 ( 1 ) m:; 


\E^L^''^ 8E^' tt'^EL^^ 8 E 2 ^ 


and V = y^x p- We have used 

^k-.w—E 




|i| 5 ^|n| 


(^2 — fi'^y 


A - 27r^ 

and q = —n. 

E 


( 22 ) 


(23) 


(24) 


When we take the large-volume limit at fixed two-pion energy, W, we expect that 2 :(s) scales as 
(and that z(0) ~ —12 ~ L^). Thus the finite volume corrections decrease as 1/L when L ^ 00 [16]. 

The terms (...;...) are the forward time contributions (ETCs) for the two cases ti ^ t 2 and t 2 = ti 
respectively. These are the terms which are reabsorbed by the sink when the matrix element is extracted 
(for a detailed discussion see section 4.1 of ref. [4]). 

From the above equations we find 


AIT 

Vll 

ETCs 


^ (o ^ . 7 

4/2L3£; 2 ) 4y2^3 


1 - 


^ [{ELf \8E^ 
3m,^ 


-1 - 


ao = 


z{l)_ 
tt^EL 


Tm-Tr 

167r/2 


1 - 


8 e ;2 


i-k 


( rzji) 


1+1- 




8E2 J \ipEL^’ 6pEL^ 


1 -k 


167r2/2 TTIttE 

5 


24/2to^L3’ ASpmTrL^ 


Stt^ \ 


(25) 


where the limits refer to the case with the two pions at rest. All these results are in agreement with 
expectations: the energy shift in a finite volume is precisely the one predicted by the Liischer quantization 
condition [ 6 ]; the LL factor is in agreement with the general formula of ref. [1] and the ETCs term will 
be cancelled when we divide the correlation function by the square root of the tttt ^ tttt correlator of 
eq. (7). 


D Finite-Volume Scalar Correlation Function in Quenched QCD 


In this appendix we give the complete one-loop expression of the finite-volume correlation function for 
the scalar operator with two pion fields in the quenched theory: 
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where y = wY/mY and mo and a are the parameters characterizing the y' propagator [4], Evaluating 
the sums, we obtain: 
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